Heat capacity of the quantum magnet TiOCl 
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Measurements of the heat capacity C(T,H) of the one-dimensional quantum magnet TiOCl are 
presented for temperatures 2 K < T < 300 K and magnetic fields up to 5 T. Distinct anomalies 
at 91 K and 67 K signal two subsequent phase transitions. The lower of these transitions clearly 
is of first order and seems to be related to the spin degrees of freedom. The transition at 92 K 
probably involves the lattice and/or orbital moments. A detailed analysis of the data reveals that the 
entropy change AS through both transitions is surprisingly small (~ 0.17?), pointing to the existence 
strong fluctuations well into the non-ordered high-temperature phase. No significant magnetic field 
dependence was detected. 

PACS numbers: 71.30.+h, 72.80.Ga, 65.40.Ba 
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I. INTRODUCTION 

The discovery of high-Tc superconductivity in the 
cuprates and of colossal magneto-resistance in manganite 
perovskites generated an enormous interest in transition- 
metal oxides (TMOs). TMOs are characterized by an 
intimate coupling of spin, charge, orbital and lattice de- 
grees of freedom which is the origin of a number of com- 
plex and exotic ground states. The discovery of a num- 
ber of new low-dimensional spin 1/2 quantum magnets 
was another result of this renewed interest. After the 
first experimental observation of a spin-Pcierls scenario 
in the organic compound TTF-CuBDTi more than 30 
years ago, CuGeOg^ was the first inorganic spin-Peierls 
system. Another paramount example is NaV2 05^ a spin 
ladder with mixed-valent vanadium ions undergoing a 
charge ordering transition . Further interest in these low- 
dimensional quantum magnets comes from the fact that 
upon doping they may undergo metal-to-insulator tran- 
sitions and possibly reveal superconductivity. During the 
last decade TiOCl has been a distinguished candidate for 
a resonating- valence bond ground state 5 and has emerged 
as a further quantum-spin magnet with a gapped ground 
stated. 

Titanium-oxochloride has first been synthesized by 
Friedel and Guerin in 1876i£ The first detailed report on 
growth and characterization has been given by Schafer 
et al£ almost 50 years ago. TiOCl crystallizes in the 
orthorhombic FcOCl structure, consisting of Ti-0 bilay- 
ers separated by CI layers. In this structure every Ti is 
surrounded by 4 oxygen and 2 chlorine ions forming a 
distorted octahedron. Originally, based on an analysis of 
susceptibility data combined with band-structure calcu- 
lations, it has been suggested that the t2 g orbitals are or- 
bitally ordered, producing one-dimensional antiferromag- 
netic S = 1/2 chains 6 . The high-temperature magnetic 
susceptibility could well be described by a Bonner-Fisher 
type of behavior with an exchange constant J = 660 K. A 
sudden drop of the susceptibility almost to zero indicated 
the opening of a spin gap at T c % — 67 K, while a further 
anomaly has been detected at T C 2 « 95 K \ 6 Electron spin 



resonance^, infrared and Raman spectroscopjiiSiii as well 
as NMR measurements^ substantiated these results and 
established a spin gap of the order of 430 K. The ori- 
gin of the gap has been identified as due to a first order 
spin-Peierls transition at T c \, as evidenced by the obser- 
vation of Ti dimerization in temperature-dependent x-ray 
diffraction^ The nature of a second transition at T c % is 
still unclear. While a broadening of the NMR lineshape 
into a wide continuum below T c i has been attributed to 
a possible incommensurate (orbital) ordering^, no cor- 
responding superlattice reflections could be observed in 
x-ray scattering^ though there is some evidence for a 
structural symmetry lowering from the high-temperature 
phase. Close to T* w 135 K the opening of a spin pseudo 
gap as detected in NMR™ and the onset of giant phonon 
anomalies in infrared and Raman spectroscopjiii indi- 
cated the presence of strong spin fluctuations and a pro- 
nounced coupling to the lattice, respectively, for temper- 
atures up to T* w 135 K. 

Electronically, TiOCl is a Mott insulator with nomi- 
nally trivalent Ti. Electronic structure calculations using 
the LDA + U predict a Ti "id l xy ground stated how- 
ever, coupling to optical phonon modes can lead to strong 
orbital fluctuations within the ti g crystal field multiplet. 
The importance of correlation effects has been elucidated 
in LDA + DMFT^L studies. 

In this brief experimental report we provide detailed 
measurements of the heat capacity of TiOCl to further 
unravel the nature of the two phase transitions. 



II. EXPERIMENTAL DETAILS 

Single crystals of TiOCl were prepared by chemical 
vapour transport from TiCl3 and Ti02. 8 The samples 
have been characterized using x-ray diffraction, SQUID 
and ESR measurements at X-band frequencies. The crys- 
tal symmetry was found to be orthorhombic with lat- 
tice parameters of a — 0.379 nm, b = 0.338 nm and 
c = 0.803 nm. The magnetic properties were found to 
be in excellent agreement with published results^. The 



2 



TiOCI 



o 



O 




300 



T(K) 



FIG. 1: Heat capacity of TiOCI vs. temperature. The ver- 
tical bars indicate experimental uncertainties and reveal the 
scatter of C(T) in different measuring cycles with different 
temperature stimuli and using different samples. Upper in- 
set: Magnified region of the two phase transitions. The heat 
capacity measured at 5 T (stars) is compared to the heat 
capacity in zero external field (circles), both measured on 
cooling. The temperature region around the lower transi- 
tion at T c i « 66.7 K was remeasured (triangles) utilizing 
a reduced relaxation amplitude of 5T ~ 0.5 K on heating 
and subsequent cooling. For clarity these data are shifted by 
50 mJ/molK 2 Lower inset: Low-temperature heat capacity 
plotted as C/T vs. T 2 . For the results in zero magnetic field 
a fit using a Debye-derived phonon contribution is indicated 
as solid line. 



heat capacity measurements have been performed with 
a commercial PPMS from Quantum Design for tempera- 
tures 1.8 K < T < 300 K and in external magnetic fields 
up to 5 T. 



III. RESULTS AND DISCUSSION 

Fig. 1 shows the central result of this investigation, the 
heat capacity as function of temperature. The heat ca- 
pacity is characteristic for a three-dimensional solid with 
a Debye temperature of order 200 K. Superimposed we 
find two weak but distinct anomalies which are clearly 
related to the two subsequent phase transitions at T c \ — 
67 K and T c2 = 91 K. In recent NMR experiments^ 
slightly different temperatures, namely 94 K and 66 K, 
respectively, have been determined. The effects of the 
phase transitions on the specific heat are remarkably 
weak and are almost lost under the large phonon- derived 
heat-capacity contributions, and it is clear that the en- 
tropy covered by the two anomalies is comparably small. 

The upper inset in Fig. 1 shows the heat capacity 
(given as C/T vs. T) in the temperature region of the 



phase transitions for both zero field and an external mag- 
netic field of 5 T measured on cooling together with a 
special heating/cooling run in the vicinity of T c \ in zero- 
field, which will be discussed later. 

Concerning the magnetic field dependence we observe 
no effects on the heat capacity within experimental un- 
certainties. The moderate magnetic fields used here nei- 
ther shift the transition temperatures, nor do they seem 
to affect the entropies involved in the phase transitions. 
On the other hand, the small energy of our fields of order 
50 kOcx/^s (« 0.06 meV) is negligibly small compared 
to the intrinsic magnetic energy scale of TiOCI as esti- 
mated from ksTci^ (5.8 and 7.9 meV) or the exchange 
constant J (57 meV). 

In order to analyze the nature of the detected anoma- 
lies in terms of first or second order phase transitions, we 
performed heating/cooling cycles across the phase tran- 
sitions. As the experimental setup utilizes a relaxation 
method, each data point is related to the average over a 
temperature interval above the initially stabilized tem- 
perature rather than to an exact temperature^! This 
means, that for the cooling sequence the temperature 
is decreased between the data acquisition and increased 
during the acquisition process itself. Thus, in the case 
of hysteretic (i.e. first order) behavior, the actual tran- 
sition, together with the corresponding release of latent 
heat, only is fully captured only in the heating branch of 
the measuring sequence. Around T c \ a significant differ- 
ence between heating and cooling can be detected, point- 
ing towards a first order transition (shifted curve in the 
inset of Fig. 1). No such feature could be found for the 
upper transition at T c i- 

The lower inset in Fig. 1 shows the low-temperature 
heat capacity plotted as C/T vs. T 2 to search for possible 
spin contributions. In zero external magnetic field, C(T) 
is fully determined by phonon contributions and from a 
linear fit (solid line) a Debye temperature of 9^, = 210 K 
can be determined. The low-temperature heat capacity 
becomes slightly enhanced in magnetic fields of 5 T which 
can be explained by the increasing importance of nuclear 
hyperfine contributions or the influence of paramagnetic 
defects at the lowest temperatures. For comparison, the 
low-temperature specific heat of the canonical organic 
spin-Peierls compound (TMTTF) 2 PF 6 1S displays besides 
the nuclear hyperfine and T 3 phonon contributions an 
additional quasi-linear (~ T 12 ) term which has been as- 
cribed to low-energy excitations. A similar contribution 
certainly can be excluded for TiOCI. 

In the following analysis we will decompose the heat 
capacity into phonon and spin (S = 1/2) contributions, 
with the former dominating the latter for all tempera- 
tures. Nevertheless, we are able to perform a straight- 
forward analysis of the thermal properties of TiOCI and 
we will show that we arrive at a rather consistent descrip- 
tion despite the large uncertainties mentioned above. 

Fig. 2 shows the heat capacity for the complete tem- 
perature regime investigated, plotted as C/T vs. T. As- 
suming a S = 1/2 spin chain the contribution of the spin 
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chain for T > 67 K is modelled by a Bonner-Fisher type 
of behavior— taking into account an exchange coupling 
J = 660 K. For the low-temperature regime, i.e. below 
the spin-Pcicrls transition at T c \ — 67 K, the expected 
exponential decrease of the spin part of the heat capacity 
is simulated by a BCS-like mean-field treatment with 
a temperature dependent gap, 2Ao = 3. 5kg Tci; and, at 
T c i, a heat capacity jump of 1.43 times the contribution 
of the spin-chain, as it was proposed for canonical spin- 
Peierls systemsi Using this approach we can calculate 
the spin-derived heat capacity without any free parame- 
ter. The result (for representation multiplied by a factor 
of 3) is shown as dash-dotted line in Fig. 2. We are 
aware that this can only be a very rough analysis, specif- 
ically having in mind that the actual gap in TiOCl seems 
to be much larger and has recently been determined as 
2A = 10 - 15k B T c i, c2 by NMR experiment^ and that 
the transition actually is of first order. The disregard of 
these experimental facts gives an underestimation of the 
jump height of the magnetic specific heat at the spin- 
Peierls transition and a too broad decay below, leading 
indeed to observable differences between the measured 
data and our calculation closely below T c \. However, the 
overall entropy balance is not affected by these shortcom- 
ings of our simple model. The total magnetic entropy at 
high temperatures sums up to Rln 2 as expected for a 
spin 1/2 system. 

The phonon system has been fitted assuming one 
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FIG. 2: Heat capacity of TiOCl plotted as C/T vs. T. The 
5=1/2 contribution including a mean-field spin-Peierls tran- 
sition as described in the text is indicated as dash-dotted line. 
For representation purposes the calculated values have been 
multiplied by a factor of 3. The fit to the phonon-derived 
specific heat (see text) is indicated as dashed line. The sum 
of both contribution is drawn as solid line and represents the 
best fit to the experimental results. The inset shows the dif- 
ference between calculated and measured heat capacities in a 
limited temperature regime. 



Debye- and two Einstein-type contributions, yielding a 
total of 5 free parameters, namely, the mean Debye (©d) 
and two Einstein temperatures (6e 12 ), and the ratio of 
Debye to Einstein modes Rjj /£ , as well as the ratio be- 
tween the Einstein modes Re 1 /e 2 - The number of de- 
grees of freedom (Nf), which of course should be nine per 
formula unit of TiOCl, was kept fixed. The experimen- 
tally determined heat-capacity values for 2 K < T < 65 K 
and 130 K< T < 300 K have been used for the fit of the 
heat capacity. For the fitting procedure we included the 
spin contributions which have been calculated parameter- 
free as outlined above. The total heat capacity, spin and 
phonon contributions are shown as solid line in Fig. 2. 
The lattice derived heat capacity is indicated as dashed 
line. Despite this oversimplified model we arrive at an 
astonishingly good description of the heat capacity over 
the complete temperature range. The parameters as de- 
termined by the best fit seem to be realistic: The char- 
acteristic temperatures, 9d = 188 K, ©^ = 352 K, 
and @e 2 = 614 K, with the ratios Rd/e — 0.30 and 
Re 1 /e 2 = 0.51 determining the relative weight of Debye 
and Einstein modes, are reliable keeping in mind that 
the optical phonon modes in TiOCl range from approxi- 
mately 100 K to 700 K with dominant modes close to our 
values for © Bl 2 piSLii The small discrepancy of the De- 
bye temperature obtained from this over-all fit to C(T), 
compared to the value derived from the low-temperature 
heat capacity only, can be explained by the presence of 
relatively low lying optical modes. 

Comparing the model heat capacity to the measured 
one in Fig. 2 we find very good agreement for the smooth 
temperature evolution for low (T < T c i) and high (T > 
T c2 ) temperatures. As clearly seen, the main deviations 
arise in the temperature region between the two phase 
transitions. The good agreement concerning the smooth 
part of the heat capacity may seem surprising in the light 
of our simplified model. On the other hand, the phonon 
part of our C(T)-model is reasonably close to a realistic 
description involving the true phonon spectrum and the 
spin part of the heat capacity bears only a small weight. 
Thus, we believe that our model yields a reliable esti- 
mate of the regular part of the specific heat, allowing us 
to separate out the anomalous contribution due to the 
phase transitions at T cl and T c2 . The inset of Fig. 2 
shows this extra heat capacity in a limited temperature 
range from 50 to 150 K. The integral over this region 
gives an estimate of the entropy AS* being released when 
going through both transitions. AS is of the order of 
0.12 R (±0.02 R) vanishingly small compared to Rln(2) 
expected for a spin S — 1/2 system or an orbital doubly 
degenerate state. The narrow peak at T c \ corresponds 
to the latent heat released at the phase transition where 
according to magnetic measurements spin dimerisation 
appears. The larger fraction of the entropy is covered 
under the phase transition at T c2 extending up to 135 K, 
exactly the temperature where in NMR experiments the 
opening of a pseudo gap has been detected— Assum- 
ing that the phonon and spin derived heat capacities are 
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correctly described, we are left with the orbital degrees 
of freedom only. However, if the orbitals undergo an or- 
dering transition at T C 2 the entropy again is much too 
low and orbital fluctuations must extend to much higher 
temperatures, even far above T*, the characteristic tem- 
perature of the opening of a pseudo gap. In this sense 
our results are compatible with the strong fluctuation 
effects observed in NMRiSi and in Raman scatteringii 
and with the results of density-functional calculations 
suggesting that TiOCl may be subject to strong orbital 
fluctuations^. This indicates that the phase above T c \ is 
characterized by the presence of pronounced fluctuations 
retaining some sort of short range order and thereby de- 
laying the release of the full entropy with temperature, as 
typical for one-dimensional systems. Unfortunately, our 
specific heat data do not allow us to decide on the pre- 
cise microscopic nature of the fluctuations (spin-Peierls 
or orbital ordering). 

IV. CONCLUSIONS 

Summarizing, we have measured the heat capacity of 
TiOCl for temperatures 2K<T< 300 K and magnetic 
fields up to 5 T. We observed two anomalies at T c \ = 91 K 
and T C 2 = 67 K corresponding to two subsequent phase 
transitions. The phase transition at 67 K, where strong 



spin-dimerization sets in, reveals significant hysteresis ef- 
fects and hence is of first order. Most probably the orbital 
degrees of freedom are involved in the phase transition 
at T C 2, which shows no thermal hysteresis. The entropy 
below this anomaly is vanishingly small indicating that 
the largest fraction of entropy is released at considerable 
higher temperatures. The total heat capacity can sat- 
isfactorily be explained describing the phonon spectrum 
with Debye and Einstein modes and the magnetic excita- 
tions by a Bonner-Fisher type heat capacity for T > T c \ 
and the opening of a spin gap below. At low temperatures 
and zero external fields the heat capacity can satisfacto- 
rily be described taking a phonon-derived T 3 term into 
account only. No indication of contribution of low-energy 
excitations can be detected. 
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